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ABSTRACT 

Dense, star-forming, cores of molecular clouds are observed to be significantly 
magnetized. A realistic magnetic field of moderate strength has been shown to 
suppress, through catastrophic magnetic braking, the formation of a rotationally 
supported disk during the protostellar accretion phase of low-mass star formation 
in the ideal MHD limit. We address, through 2D (axisymmetric) simulations, the 
question of whether realistic levels of nonideal effects, computed with a simplified 
chemical network including dust grains, can weaken the magnetic braking enough 
to enable a rotationally supported disk to form. We find that ambipolar diffusion, 
the dominant nonideal MHD effect over most of the density range relevant to 
disk formation, does not enable disk formation, at least in 2D. The reason is 
that ambipolar diffusion allows the magnetic flux that would be dragged into 
the central stellar object in the ideal MHD limit to pile up instead in a small 
circumstellar region, where the magnetic field strength (and thus the braking 
efficiency) is greatly enhanced. We also find that, on the scale of tens of AU 
or more, a realistic level of Ohmic dissipation does not weaken the magnetic 
braking enough for a rotationally supported disk to form, either by itself or in 
combination with ambipolar diffusion. The Hall effect, the least explored of 
these three nonideal MHD effects, can spin up the material close to the central 
object to a significant, supersonic rotation speed, even when the core is initially 
non-rotating, although the spun-up material remains too sub-Keplerian to form 
a rotationally supported disk. The problem of catastrophic magnetic braking 
that prevents disk formation in dense cores magnetized to realistic levels remains 
unresolved. Possible resolutions of this problem are discussed. 
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Introduction 



The formation and early evolution of disks is a long-standing fundamental problem in 
star formation. Early work in the field had concentrated on the simpler problem of disk 
formation fr om the collap s e of a rota ting d e nse c ore in the absence of a magnetic field, as 
reviewed in iBodenheimerl (jl995[ ) and iBossI (jl998l ). Dense star- forming cores are observed 
to be significantly magnetized, however. There is increasing theoretical evidence that disk 
formation is greatly modified, perhaps even suppressed, by a dynamically important magnetic 
field. 

The most comprehensive measurer nents of the magneti c field strength in dense cores 
of low-mass star formation come from iTroland fc Crutcherl (120081 ). They carried out an 
OH Zeeman survey of a sample of nearby dark cloud cores, probing densities of order 10^- 
lO^cm"^. The inferred m ean value for the diraensionless mass-to-fiu x ratio, relative to the 
critical value jNakano fc NakamuralEgTsl . Ishu fc Li|[l997l ). is Xios ^ 4.8 ± 0.4. It 

was obtained from the measured line-of-sight component of the magnetic field Bios without 
any geometric correction. Correc ting for geometr ic effects statistically would lower the mass- 
to-flux ratio by a factor of 2-3 ( IShu et al.l 119991 ). bringing the mean value of the intrinsic 
mass-to- flux ratio to a few, i.e., A ~ 2-3. Such values of A are naturally produced in 
the scenario of ambipolar diffusion-regulated dense core formation in strongly magnetized 
(mag netically subcritical) clouds, with or without the assistance of supersonic turbulence 



Lizano fc Shu 1989: Basu fc Mouschovias 1994: Nakamura fc Li 2005: Kudoh fc Basu 



201ll ). It may also be consistent with the scenario of dense core formation from turbulence 
compression in more weakly magnetized background clouds, since the core is expected to 
be more strongly magnetized relative t o its mass (i.e., lower A) than the cloud as a whole 
( iTilley fc Pudritzll2005l : iDib et al.ll2007l ). Well ordered magnetic fleld s are also inferred from 



polar ization maps of du st continuum emiss ion, both on the core scale (jWard-Thompson et al. 



2000l ) and smaller (e.g.. iGirart et al.ll2006l ). We therefore expect the dense cores to be rather 
strongly magnetized based on both observational data and core formation theories. 

A moderately strong magnetic fleld can suppres s disk formation in the ideal MHD 
limit. This was flrst demonstrated in lAUen et al.l ( 120031 ). who carried out 2D (axisymmetric) 
simulations of the collapse of rotating cores magnetized to a level of A < 10. The basic reason 
for disk suppression is that, in the ideal MHD limit, flux freezing allows the infalling material 
to drag a flnite amount of magnetic flux into the central object, creating a split magnetic 
monopole whose (poloidal) fleld strength increases rapidly with decreasing radius (as r~^, 
Galli et al.ll2006l : see their Fig. 1). The increased fleld strength close to the central object, 
coupled with a long magnetic lever arm from severe equatorial pinching of (poloidal) fleld 
lines, is responsible for the catastrophic magnetic braking that suppresses the formation of a 
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rot ationally sup p orted disk. This magnetic braking; catast rophe was confirmed numerically 
by iMellon fc Lil (I20081) and iHennebelle fc Fromand ( l2008l ) using, respectively, 2D and 3D 



simulations (see also Price fc Bate 



aligned (see, however, iMachida et al. 



2007D. a t least when the magnetic and rotation axes are 



2010 and discussion in 



When the magnetic and rotation axes are misaligned, iHennebelle fc Ciardil ( l2009l ) found 
in their AMR MHD simulations that magnetic braking efficiency is reduced relative to the 
aligned case. A potential concern is that the misalignment increases the flow complexity, 
which may enhance the nur nerical magnetic diffusio n that is considerable on small scales (see 
a nice discussion in §5.3 of IHennebelle et al.ll201ll ). A more obvious possibility for avoiding 
the magnetic braking catastrophe is to relax the ideal MH D approximation. Since dense cores 
are known to be lightly ionized f lBergin fc Tafallal 120071 ). no n-ideal MHD effects (including 



ambipolar diffusion. Hall effect and Ohmic dissipation, e.g., iNakano et al.ll2002l ) are to be 
expected. The fi rst non-ideal MHD effect considered in this context was Ohmic dissipation. 
Shu et al.l (120061 ) suggested that in order for Ohmic dissipation to weaken the field strength 
over a large enough region so that a large-scale rotationally supported disk of tens of AUs 
or more can potentially form, the resistivity must be at least one order of magnitude above 
the cl assical (microscopic) value. This suggestion was confirmed by iKrasnopolsky et al. 
(120 lOl ). who found that large, 100 AU scale disks can indeed form, as long as t he resistivity is 



enhanced by a large enough factor, to a value of order 10^^ cm^ s 



Machida et al. 



or more 

(I2OIOI ) carried out core collapse calculations including a distri bution of resistivity w ith density 
and temperature (from a fit to the resistivities computed in lNakano et al.ll2002l ) and found 
that, even with just the classical resistivity, a small rotati onally supported dis k can form at 
the beginning of the protostellar accretion phase (see also lDapp fc Basull2010l) and grow to 
l arger , 100- AU scales at later tiines. Pa rt of the apparent discrepancy between iMachida et al. 
(120 lOl ) and IKrasnopolsky et al.l (120 lOl ) may be due to different simulation setup. Another 
difference may be in the level of numerical magnetic diffusivity (see discussion in Further 
investigation is needed to clarify the situation in the limiting case that includes only Ohmic 
dissipation. 



Ohmic di ssipation is important at h igh densities (above ~ 10 cm , iNakano et al. 



2OO2I : see also iKunz fc MouschoviasI l2010l ). Before reaching such densities, dense core ma- 
terial must evolve through lower densities, where the Hall effect and especially ambipolar 
diffusion dominate over Ohmic dissipation. The effect of ambipolar d iffusio n on disk for- 
mation was investigated s emi-anal yticallv by IKrasnopolsky fc Konigll (|2002|). and riumer i- 



cally bv lMellon fc Lil ( l2009l : see also lDuffin fc Pudritzll2009l and lHosking &: Whitworthll2004h . 



Compared to the ideal MHD case, a new ingredient is the ambipolar diffusion-induced 
accretion shock, driven by the magnetic f lux left behind by th e material that has gone 



into the central object (iLi fc McKeelll996l : ICiolek fc Konigllll998l : IContopoulos et al.lll998 
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Tassis fc MouschoviasI 120071 ). iKrasnopolsky fc Konigll (j2002[ ) demonstrated using a ID semi- 
analytic model that the strong magnetic field piled up inside t he ambipolar diffusi on (AD) 
shock can in principle brake the post-shock material efficiently. iMellon fc Lil ( 120091 ) showed 
through 2D (axisymmetric) simulations that this is indeed the case. They started thei r 
calculations from a self-similar rotating, magnetized isothermal toroid (lAllen et al.ll2003l ). 
A power-law dependence on the neutral density is assumed for the ion density, so that the 
subsequent collapse remains self-similar, even in the presence of ambipolar diffusion. The 
self-similarity provides a powerful check on the validity of the numerically obtained collapse 
solutions. It imposes, however, strong restrictions on both the initial core properties and 
charge densities. 

One objective of the present paper is to investigate the ro le of ambipo l ar diff usion 
in disk formation without the restrictive simplifications made in iMellon fc Lil (j2009[ ). We 
do this by following both the pre-stellar evolution of the rotating, magnetized dense core 
and the protostellar mass accretion phase after a central stellar object has formed and by 
computing the charg e densities self-cq n sistently using a si mplified chemical network that 



includes dust grains (INishi et al. 



1991 



Nakano et al. 



20021 ): they can affect the magneti c 



diffusivities greatly (e.g.. IWardle fc Nglll999l ). Another objective is to extend IMellon &: Li 



(120091 ) calculations by including both Ohmic dissipation and Hall effect in addition to am- 
bipolar d iffusion. The Hall effe ct was e xplored previously in the context of disk-driven 
outflows (IWardle fc Konigll Il993l : see also iKonigl et al.l |2010[ ) and accretion disk dynamics 



(ISano fc Stone 



2OO2I): it is only starting t o be explored in t he context of core collapse and 



disk formation (IKrasnopolsky et al.l 1201 It IC. Braidingll201l[ ). 



We find that, on scales greater than 10^^ cm (or 6.7 AU) that we are able to resolve 
in our non-ideal MHD calculations, no rotationally supported disks form in dense cores 
with a moderately strong magnetic field (with A ~ several), largely because of the excessive 
braking due to the magnetic field trapped interior to the AD shock. On these scales, Ohmic 
dissipation affects the flow dynamics relatively little. The Hall effect, on the other hand, 
can torque up the spun-down post AD-shock material to signiflcant, supersonic, rotation 
speeds. The rotation speeds remain well below Keplerian and rotationally supported disks 
are not formed. Our non-ideal MHD calculations re-enforce the idea that disk formation 
is difficult in the presence of a moderately strong magnetic fleld, at least in 2D (assuming 
axisymmetry) . We discuss possible ways to get around this difficulty in the discussion section 
(©• 

The rest of the paper is organized as follows. In §2], we describe the problem setup, 
including the initial and boundary conditions, the nonideal MHD code used in this work, and 
the computation of charge densities including dust grains. Disk suppression by a moderately 
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strong magnetic field under a range of realistic conditions is demonstrated in ^ We discuss 
potential disk formation in the case of weaker magnetic fields in §1] and the spin up of an 
initially non-rotating collapsing envelope by the Hall effect in §3 Our main results are 
summarized §71 



2. Problem Setup 



2.1. Initial and Boundary Conditions 

Low-mass pre-stellar core s in nearby star-forming regions are observed to have relatively 
simple dynamical structures f lBergin fc Tafallal 120071 ) . We idealize such cores as initially 



unifo rm spheres oi R = 10^'' cm in mass and 1 Mq in mass (see 
20081 for a similar setup). The initial core mass density is thus po = 
corresponding to a volume density for molecular hydrogen = 10^ cm~'^ 
hydrogen nuclei for each He), and a free fall time t// 
equation of state is adopted below a critical density p 



Hennebelle fc Fromang 



a = 0.2 km s ^) and P oc p'^/^ at densities above (e.g., 
of thermal to gravitational binding energy is a 



4.77 X 10-^9 gcm-3, 
assuming 10 
3 X 10^^ s. A simple isothermal 
10~^'^gcm~^ (with a s ound speed 



Masunaga &: Inutsukall2000l ) . The ratio 



2.5Ra^/{GM) = 0.75. On this uniform 
sphere, we impose a uniform magnetic field Bq at the beginning of the simulation. We choose 
a fiducial value for Bq of 10~^ in the Lorentz-Heaviside units that are convenient for the Zeus 
family of codes. It corresponds to Bq = y/Anx 10~^ G = 35.4 pG in Gaussian CGS units. The 
ratio of magnetic to gravitational binding energy is 7 = 0.13. Another way to characterize the 
strength of the magnetic field is through the dimensionless mass-to-flux ratio A, in units of the 
critical value (27rG'^/^)~^. For the core as a whole, the mass-to-fiux ratio is A = 2.92. On the 
central flux tube that passes through the origin, the mass-to-fiux ratio is Ac = 4.38, higher 
than the global value by 50%. These fiducial values of mass-to-fiux ratios are somewhat 
larger than the values of A ~ 2-3 obtained in models of ambipola r-diffusion driven core 



formation out of ma gnetically su bcritical background clouds, with (iNakamura &: Lil 12005 



Kudoh &: Basull201ll ) or without flLizano fc Shijll989l : iBasu fc MouschoviasI 119941) turbulent 
compr ession. They are consistent with the mean value of A inferred by iTroland fc Crutcher 
(120081 ) after correcting for geometric effects (see §1). 



One may argue that the adopted fiducial field strength of Bq = 35.4 pG at a number 
density = 10^ cm~^ is too high, because such a field strength is rarely measured in 
nearby regions of low-mass star formation (such as the Taurus clouds) directly using OH. 
However, the OH ernission is dominated by relatively low density envelopes of dense cores 
(ITroland &: Crutcherll2008l ). If the core material was condensed out of a more diffuse gas, its 
field strength would be lower before condensation. For example, at densities of order ~ 
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3.5 3 (^probed by typical OH observations; ICrutcher et al.ll201Cll ). the pre-condensation 



10 

field strength of our model core would be only ~ 3.5 fiG (if the condensation occurs more or 
less isotropically under flux-freezing condition, as it should be since the core is substantially 
magnetically supercritical). It is lower than the me dium value of ~ 6 /iG inferred for the 



atomic cold neutral medium (IHeiles fc Trolandll2005l ). If anything, our adopted fiducial field 
strength may be on the low side. Nevertheless, we will consider some smaller values for Bq 
as well, in view of the possibility that t here may be a broad distribution of field strengths 
at a given density ( ICrutcher et al.l 120101 ) and the fact that the majority of stars are formed 
in clusters where the magnetic field is less well observed. 

One may also be concerned that our adopted initial d ensity distribution is u niform, 
whereas the observed pre-stellar cores are centrally condensed ( iBergin fc Tafallall2007l ) . How- 
ever, the density distribution changes with time, evolving through a series of centrally con- 
densed configurations with different degrees of central-to-edge contrast, as seen in Fig. [T]in 
the absence of magnetic field and rotation. For example, the central density n^a ~ 10^ cm~^ 
at time 2.5 x 10^^ s is close to that inferred for the well studied pre-stellar core L1544 
( IWard-Thompson et al.lll999l ). although the maximum infall speed at this time is somewhat 
highe r than the observationally inferred value, which is roughly lO^cms"^ ( iTafalla et al. 
19981 ). The infall speed is reduced, however, by both rotation and especially magnetic fields 
that are included in the majority of our calculations. Indeed, the infall is approximately 
half-sonic for our reference model (see Table 1) when the central density ~ lO^cm"^, in 
agreement with observations. 

An alternative approach would be to start with a centrally condensed static configuration 
that mimics the obser ved pre-stellar densi ty profile, with a uniform magnetic field imposed 
at the beginning (e.g.. iMachida et al.ll2007l ). Implicit in this approach is the assumption that 
there is strong evolution in the density profile of the core but no corresponding evolution 
in its magnetic field distribution. For our purpose of studying the efficiency of magnetic 
braking, it has the drawback that the central flux tubes are loaded with much more mass 
than those in the outer part. It would lead to a weaker magnetic braking in the early phase 
of protostellar collapse that, in our view, is hard to justify physically. 

The rotating profile of dense cores is not well constrained by observations. For simplicity, 
we adopt for our initially uniform core a solid body rotation, with a fiducial angular speed 
Q = 10~^^s~^. It corresponds to a ratio of rotational to gravitational binding energy /3 = 
0.025, which is typical of t he values inferred for NH3 cores based on the velocity gradient 
observed across the cores (iGoodman et al.l 119931 ): the inferred /3 has a considerable range, 
which motivates us to consider other values of /3 as well. 



As in iMellon fc Lil (120081 12009[ ). we adopt a spherical polar coordinate system {r,6, 



10" 



1 0' r 



0) 



-2x10' 




10'" 10' 
Radius (cm) 



10 



-3x10* 



Radius (cm) 



Fig. 1. — Distributions of the density (left panel) and infall velocity (right) at three times 
2 X 10^^ (dotted line), 2.5 x 10^^ (solid) and 3 x 10^^ s (dashed) for the collapse of a non- 
rotating, non-magnetize dense core. 
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for our axisymmetric simulation. The inner boundary is set at (or 6.7 AU) and the 

outer boundary at 10^'' cm. Even though our inner radius is relatively large compared with 
other collapse studies, such as Machida et al. (2010) who set the sink particle size to 1 AU, 
our smallest cell size is only 0.2 AU, which ensures that the flow dynamics near the inner 
boundary is well resolved. The high resolution is particularly important for minimizing the 
numerical magnetic diffusion that may affect the trapping of magnetic flux at small radii, 
which lies at the heart of efficient magnetic braking and disk suppression. The standard 
outflow boundary conditions are enforced at both the inner and outer boundaries. Matter 
that crosses the inner boundary is collected at the origin. The central point mass interacts 
with the matter in the computation domain through gravity. We use a computational grid 
of 120 X 90 that is non-uniform in the r-direction, with a spacing Ar = 0.2 AU next to the 
inner boundary. The spacing increases outward by a constant factor ~ 1.0647 from one cell 
to the next. The grid is uniform in the ^-direction. 



2.2. Induction Equation and Non-Ideal MHD Code 



At the heart of our non-ideal MHD core collapse problem lies the induction equation 



9B f B 

— = Vx(vxB)-Vx[r/o(VxB)]-Vx|7/H[(Vx B) x -] 



-Vx<;r/^|x[(VxB)x|: 



where the Ohmic, Hall and ambipolar diffusivities are related to the electric con ductivity par- 



allel to the field line a\\, Pedersen and Hall conductivities ap and an through (INakano et al. 
2OO2I ) 

(2) 



Allan 



Vh 



Va = -- 



47r \aj, + ajj 



An al + ajj' 

These diffusivities will be discussed further in the next subsection. Here, we describe briefly 
our numerical treatment of the three non-ideal MHD terms. 



Our code, dubbed "ZeusTW," was derived from the ideal MHD code Zeus3D fIClarke et al 



1994f). It treat s the O hmic ter m in the induction equati on (1) using an algorithm based on 
Fleming et al.l (120001 see also iKrasnopolskv et al.l 120101). A mbip olar diffusion was t reated 
using the explicit method described in lMac Low et al.l ( 19951 ). as in lMellon fc Lil ( 20091 ). The 
magnetic field is evolved using a velocity that is the sum of the bulk neutral velocity v and 
a "drift" velocity defined as 



VA_ 

52 



(V X B) X B. 



(3) 



In the widely discussed limit where ambipolar diffusion is the dominant non-ideal effect and 
ions are well tied to the magnetic field, the velocity A\ad can be interpreted as the ion- 
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neutral drift velocity. In the more general case that we are studying, this is not necessarily 
the case. Nevertheless, we will define an "effective ion velocity" 



V.eff = Av^£, + V, 



(4) 



to make contact with previous work. It provides a useful measure of the effect of ambipolar 
diffusion even in the prese nce of other non-idea l M HD effects. L astl y, our treatment of 



the H all term was based on ISano &: Stond (120021 ) and iHubal (120031 . see iKrasnopolsky et al. 



20 111 for detail). The time steps for evolving the ambipolar diffusion and Hall terms in the 



induction equation are particularly stringent near the polar axes, where the magnetic field is 
strong but the density is often low, due to either gravitational collapse along the field lines 
or outfiow. In some cases, a fioor is imposed on the time step, by decreasing the ambipolar 
or Hall diffusivity in a small volume or increasing the density in an evacuated region to limit 
the Alfven speed. We have verified that the fioor has little effect on the flow dynamics. 



2.3. Charge Densities and Magnetic Diffusivities 



The magnetic diffusivities depend on the densities of charged particles, including molec- 
ular and atomic ions, electrons, and charged d ust grains. We will follow Nakano and collabo- 



rators in computing the charge densities (e.g., iNakano et al.ll2002l : iNishi et al.lll99l[ ). using a 



simplifled chemical network and simple prescriptions for grain size distribution. The network 
includes neutral species H2, He, CO, O, O2 and other heavy metals (denoted collectively by 
"M") and charged species e~, H+, He"*", C"*", M"*", B.^, and other molecular ions (denoted 
colle ctively by "m^") , as well as neutral and positively and negatively charged dust grains 



see 



Nishi et al.l Il99ll for detail). The grain size distribution in dense cores of molecular 



clouds is relatively unconstrained observationally. For illustration, we will consider the stan- 
dard MRN distribu tion (dn/da) oc a~^'^ with the grain size a between 5nm < a < 250 nm 
( iMathis et al.l Il977l ) and a grain mass that is 1% of the total. It is likely, however, for 
the grains to grow substantially in dense cores of molecular clouds; the MRN distribution, 
which is more appropriate for diffuse clouds, may be regarded as the starting point for grain 
growth in dense cores. Dire ct evidence for grai n growth comes from the Spitzer detection 
of the so-called "coreshine" ( jPagani et al.ll2010l ). which indicates that at least some grains 
have grown to micron-size. To illustrate the effects of grain growth, we will also consider an 
opposite limit where the grains have a single, large, size a = 1 fim (denoted LG distribution). 
The MRN and LG distributions should bracket the real situation where some grain growth 
is expected. 

Fig. |2] plots the fractional abundances (relative to the number density of hydrogen 
nuclei ri//) of the charged atomic and molecular species and dust grains as a function of 
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uh (related to mass density through p = 2.34 x 10~'^'^nH), with the cosmic rate ionization 
rate ( normahzed to the standard value 10~^^s~^; the normalized rate is denoted by = 
C/(10~^^ s~^). In the (simpler) large grain (LG) case, the dominant charges are metal ions 
(M"*") and electrons (e~) over the whole density range of interest to us. In the MRN case, the 
charge densities are lower compared to the LG case, because of a large amount of small grains, 
which provide a large total surface area for ions and electrons to recombine on. The metal 
ions and electrons remain the most abundant charges at low densities (below ~ 10^ cm^'^). At 
higher densities, negatively and positively charged (small) grains become more dominant. 
The MRN case w as computed with 10 size bins equally spaced logarithmically, following 



Nishi et al.l (jl99ll ). When 20 bins are used, the charge densities change by less than 3% over 



the density range that spans more than 10 orders of magnitude. 

To obtain the magnetic diffusivities in the induction equation, we need not only charge 
densities, but also magnetic field strength. The field strength will be computed self-consistently 
in our MHD simulations. The rough magnitude of the diffusivities can be illustrated using 
e.g. the field strength-density relation 

B = 1.43 X 10"^n^/^ (5) 



assumed in lNakano et al.l (120021 ). relevant for a magnetically supported sheet in hydrostatic 
equilibrium along the field lines. The computed diffusivities are shown in Fig. [3J As expected, 
the ambipolar diffusivity dominates at relatively low densities for both MRN and LG cases. 
For the MRN case, which has a large amount of small grains, the Hall diffusivity becomes 
comparable to the ambipolar diffusivity over a wide range of density (~ 10^ - lO^^cm"^). 
It has a negative value, dominated by the contribution from negatively charged grains, 
although the contribution from positively charged grains approaches (from below) that from 
the negative grains at densities greater than ~ 10^°cm~'^. At densities greater than ~ 
10^^ cm"'^ (which is generally not reached in our computational domain), Ohmic diffusivity 
dominates. The large Ohmic diffusivity r]o is the reason for the ambipolar diffusivity tja to 
become negative above a density ~ 10^^ cm~^, since it contributes negatively to rjA (see eq. 
[2]). In the opposite extreme of large grain case (LG), where all small grains are removed, 
the Hall diffusivity is positive, dominated by the contribution from metal ions. It exceeds the 
ambipolar diffusivity above a density uh ~ 10^^ cm~^. In this case, the level of ionization is 
high enough that the Ohmic diffusivity remains unimportant throughout the plotted density 
range (up to ~ lO^^cm"^). Overall, the elimination of small grains in the LG case makes 
the magnetic field better coupled to the bulk neutral material compared to the MRN case, 
as expected. 
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Fig. 2. — The number densities of charges for the standard MRN grain size distribution (left 
panel) and large, 1 /im grains (right). 




Fig. 3. — The absolute values of the computed magnetic diffusivities for the MRN (left 
panel) and large grain (right panel) cases, for the illustrative magnetic field-density relation 
given in eq. (5). For the MRN case, the Hall diffusivity rjH is negative everywhere and the 
ambipolar diffusivity tja changes from positive to negative as the density uh increases above 
~ lO^^cni-s. 



-VI- 



3. Disk Suppression by a Moderately Strong Magnetic Field 

In this section, we will consider a magnetic field of moderate strength i?o = 35.4 /iG, cor- 
responding to a dimensionless mass-to-flux ratio of A = 2.92 for the dense (n^a = lO^cm"^) 
core as a whole and 4.38 for the central flux tube. As we argued in this field strength 
is likely on the conservative side. Since the ambipolar diffusivity dominates over the Ohmic 
and Hall diffusivities in most of the density range encountered in our calculations (see Fig. 
|3]), we will first concentrate on its effect on core co llapse and disk forrn ation in §3. H and §3.21 



Similar AD-only calculations were performed by iMellon fc Lil (120091 ). except that we now 
include the prestellar phase of core evolution leading up to the formation of a central object 
(in addition to the protostellar phase that they studied) and a more detailed calculation of 
the charge densities, including charged grains. The additional effects of the Ohmic and Hall 
diffusivities are considered in §3.31 and §3.4[ respectively. 



3.1. Reference Model with Ambipolar Diffusion Only 

For the reference model, we adopt the fiducial values for the initial field strength 
-Bq = 35.4 /iG, cosmic ray ionization rate C = 10~^^s~^, and initial core rotation rate 
f2o = 10~^^s~^, as well as an MRN grain size distribution. Since grain growth in dense 
cores tends to make the magnetic field better coupled to the neutral matter (see Fig. [3]) 
and magnetic braking more efficient, our adoption of the MRN distribution is also on the 
conservative side. The parameters for this (Model REF) and other models are listed in Table 
1. The results are shown in Figs. lU Eland [Hi 

Fig. H] displays the density distribution and velocity field on the meridian and equatorial 
planes for the inner part (10^-AU scale) of the computation domain, at a representative time 
t = 6 X 10^^ s (or twice the initial free-fall time), when 0.57 Mq (or 57% of the initial core 
mass) has fallen to the center. From the left panel, it is clear that the density distribution 
on the meridian plane is highly flattened, especially at high densities. The dense, flattened, 
equatorial structure is not a rotationally supported disk, however. Direct evidence against 
such a disk comes from the right panel, which shows a transition from an outer region of 
rapid rotating-infall to an inner region that is neither collapsing nor rotating rapidly. The 
transition is shown more quantitatively in the left panel of Fig. [5l where the infall and 
rotation speeds on the equator are plotted. The equatorial infall is initially slowed down 
near a relatively large radius r = 5 x 10^^ cm. It corresponds to the edge of the magnetic 
bubble inflated by magnetic braking (not shown in Fig. Hj), where a magnetic barrier forces 
the collapsing mat erial over a large solid angle into a narrow equatorial channel (see Fig. 2 of 



Mellon fc Lill2008l and associated discussion). Upon passing through the barrier, the material 
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Table 1. Model Parameters 



Mo del 


Grain 


Bq (/iG) 


C(10 's ) 


iZo(10 S ) 


RoU I 


REF 


MRN 


35.4 


1 


1 


no 


LG 


LG 


35.4 


1 


1 


no 


LoCR 


MRN 


35.4 


0.1 


1 


no 


HiCR 


MRN 


35.4 


3 


1 


no 


LoROT 


MRN 


35.4 


1 


0.5 


no 


HiROT 


MRN 


35.4 


1 


2 


no 


Xi Jl;_r 


MRN 


35.4 


1 


1 


no 




MRN 


35.4 


1 


1 


no 




MRN 


35.4 


1 


1 


no 




MRN 


35.4 


1 


1 


no 


WREF 


MRN 


10.6 


1 


1 


no 


WLG 


LG 


10.6 


1 


1 


no 


WHiCR 


MRN 


10.6 


3 


1 


no 


WLoROT 


MRN 


10.6 


1 


0.5 


no 


WLoCR 


MRN 


10.6 


0.5 


1 


yes? 


WHiROT 


MRN 


10.6 


1 


2 


yes? 


VWREF 


MRN 


3.54 


1 


1 


yes? 


NoROT:^^^ 


MRN 


35.4 


1 





no 


NoROT+^o 


MRN 


35.4 


1 





no 



Note. — a) The subscript "A," "H" and "O" denote, respectively, "ambipolar diffusion," 
"Hall effect," and "Ohmic dissipation" included in a model; models without a subscript 
include only ambipolar diffusion. The superscript "-" ( "+" ) denotes an initial magnetic field 
anti-parallel (parallel) to the rotation axis, b) RSD stands for "rotationally supported disk." 
In Model WREF, a small RSD is formed temporarily at early times; it is subsequently 
suppressed by magnetic braking. For Models WLoCR, WHiROT and VWREF, an RSD 
forms early but whether it would survive to later times remains unclear because of numerical 
difficulty. 
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resumes rapid radial infall, spinning up as it collapses, until a second barrier is encountered 
at r ~ 8 X 10^^ cm. This second barrier is induced by ambipolar diffusion, which enables the 
magnetic field lines to decouple from matter and pile up outside the central object as the 
matter accretes onto the protostar (see discussion in §!]). The piled- up magnetic field drives 
a C-shock, which slows down the accretion flow to a subsonic speed. The slowdown of infall 
leads to a high density in the postshock region, which is clearly visible in the density maps 
(Fig. Hj). Upon passing through the shock, the equatorial material re-accelerates towards the 
central object, reaching a highly supersonic infall speed at small radii. The supersonic infall 
clearly indicates that a rotationally supported disk (RSD hereafter for short) is not present. 

The lack of an RSD is even more obvious from the rotation speed on the equator plotted 
in Fig. O As the infalling material enters the AD shock, its rotation speed drops quickly. 
Over most of the postshock region, the rotation speed is nearly zero, indicating an efficient 
braking of the material that accretes onto the central object through the equatorial region. 

Why is it that the rotation of the infalling material is braked almost com pletely? The 



answer lies in the strong ma.gnetic field trapped interior t o the AD shock (ILi fc McKee 



19961 : IContopoulos et al.lll998l : iKrasnopolsky fc Koniglll2002l ). The field trapping is already 



evident in the left panel of Fig. |H which shows a pileup of poloidal magnetic field lines 
at small radii. It is shown more quantitatively in the right panel of Fig. [5l where the 
distribution of the magnetic pressure (S^/[87r]) with radius is plotted along the equator. 
Note the rapid increase in magnetic pressure (and thus field strength) across the AD shock 
(near ~ 8 x 10^^ cm). There is a corresponding drop in the ram pressure (pv^), indicating 
that the strong post- shock magne tic field is trapped by the ram pressure of the pre-shock 



collapsing flow (ILi fc McKeelll996l ). The further rise in the ram pressure at smaller radii is 
due to the gravity of the central object, which re-accelerates the equatorial material (that was 
temporarily slowed down to a subsonic speed behind the AD shock) to a highly supersonic 
speed. The "effective ion speed" in the post-shock region is much lower than that of the bulk 
neutral material (see the right panel of Fig. HJ. The relatively large drift velocity Av^/j, 
driven by a large outward magnetic force, is what enables the magnetic flux to accumulate 
in the postshock region in the first place: as more and more matter accretes onto the central 
object, more and more magnetic flux is left behind. Indeed, by the time shown in Figs. H] 
and about half of the total magnetic flux of the initial core is confined within the shock 
radius ~ 8 x 10^^ cm. The postshock region, which contains a small fraction (~ 4%) of the 
total mass, is so strongly magnetized that it becomes highly magnetically subcritical (with 
a local dimensionless mass-to-fiux ratio A well below unity) even though the whole core was 
significantly magnetically supercritical to begin with; the small infall speed (much below free 
fall value) and nearly vanishing rotation speed indicate that the material in the region is 
essentially magnetically supported. The strong postshock field is only part of the reason for 
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Fig. 4. — Density distribution (color map) and velocity field (white arrows) of the reference 
model (Model REF) in the meridian (left panel) and equatorial (right panel) planes, at 
a representative time t = Q x 10^^ s. The highly flattened, dense equatorial structure is 
not a rotationally supported disk, but rather a magnetically supported, nearly non-rotating 
pseudodisk. Also plotted in the left panel are poloidal field lines, with the same magnetic 
flux between adjacent lines. The color bars above the panels are for log(p), with gcm~^ and 
cm as the units for p and length. 



- 16 - 




Fig. 5. — Left panel: equatorial infall (lower solid curve) and rotation (upper solid) speed for 
the bulk neutral material in the reference model (Model REF) at t = 6 x 10^^ s. Also plotted 
for comparison are the Keplerian speed (upper dashed) and free fall speed (lower dash- 
dotted) based on the central mass, the sound speed (horizontal dashed line), and zero speed 
line. The rapid deceleration of infalling material near the radius ~ 8 x 10^^ cm corresponds 
to a C-shock induced by ambipolar diffusion. The "effective ion speeds" defined in eq. (4) are 
shown as dotted lines. Right panel: magnetic (solid line) and infall ram pressure (dashed) 
along the equator in cgs units, showing a rapid increase in field strength near the AD shock 
and a corresponding drop in ram pressure. 
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efficient braking. Another part is that the postshock field naturally bends outwards due to 
strong equatorial pinching (see the left panel of Fig. H]), which increases its lever arm and 
thus the braking efficiency. 

Even though the equatorial rotation speed at small radii remains close to zero at late 
times (as shown in Figs. |4]and|5]), it does reach a substantial, supersonic, value for a brief 
period of time during the transition from the phase of prestellar collapse to the protostellar 
accretion phase. Fig. E] displays the equatorial infall and rotation speeds during the tran- 
sition, for 6 times separated by At = 4 x 10^° s. At the earliest of the times shown, the 
material at small radii remains static, as expected during the pre-stellar phase of core evo- 
lution for the region within one thermal Jeans length of the origin. When the central Jeans 
length shrinks to inside our inner boundary, rapid protostellar accretion ensures, leading 
to a quick spin-up of the collapsing material. The rotation speed reaches a value as high 
as ~ lkms~^ (or ~ 5 times the sound speed) near the inner boundary before decreasing 
back again. The rapid, transient spin-up i s a fe ature that is not ca ptured by the self-similar 



solutions of iKrasnopolsky fc Konigll (120021 ) and lMellon &: Lil ( l2009l ). The spindown is clearly 
associated with the development of an ambipolar diffusion induced accretion shock (see the 
curves for infall speed), which strengthens as it propagates outward. The material in the 
post-shock region is so strongly braked that it rotates backwards for a short period of time, 
before settling down to the nearly non-rotating state shown in Figs. H] and [51 The strong 
magnetic braking prevents any rotationally supported disk larger than 10^^ cm (the size of 
our inner boundary) from forming at any time during our (long) simulation, which lasted 
until t = 9 X 10^^ s, when 90% of the core mass has been accreted by the central object. 



3.2. Pure AD Runs: Grain Size, Cosmic Ray Ionization and Rotation Rates 

In the reference model, the magnetic braking has clearly suppressed the formation of 
a rotationally supported disk (although a highly flattened, dense, magnetically supported, 
nearly non-rotating pseudodisk was formed). In this subsection, we explore how robust 
this result is, by varying the physical quantities that may affect disk formation in a lightly 
ionized, rotating, magnetized molecular cloud core. These include the grain size distribution 
and cosmic ray ionization rate, both of which affect the charge densities and thus the degree 
of coupling between the magnetic fleld and the bulk neutral matter, as well as the rate of 
core rotation. We discuss models where these quantities are varied over a reasonable range 
(see Table 1). 

As discussed in §2.3[ the size distribution of the grains in dense cores of molecular 
clouds is uncertain. We have adopted in our reference model the standard power-law MRN 
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Fig. 6. — Infall and rotation speeds along the equator for the reference model (Model REF) 
during the transition from the pre-stellar phase of core evolution to the protostellar accretion 
phase, between t = 4.4 x 10^^ and 4.6 x 10^^ s, showing the development of the AD shock 
(lower curves) and the associated strong braking of the post-shock material (upper curves). 
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distribution that includes a large amount of small grains. In this subsection, we will consider 
an opposite limit, Model LG in Table 1, where all grains are assumed to be 1 /im in size 
and the small grains are completely absent. Be sides grain size , the charge densities are 



also affected by the cosmic ray ionization rate (. iPadovani et al.l ( 120091 ) compiled from the 
literature the values of ( inferred for clouds of a wide range of column densities, from diffuse 
clouds to massive protostellar envelopes (see their Fig. 15). Most of the inferred values are 
above our reference value ( = 10^^^ s^^, although there are a few exceptions. We will consider 
a model with ( = lO^^'^s^^ (Model LoCR in Table 1), which is a lower limit to the inferred 
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values. For illustration, we will also consi der a case with a higher rate C = 3 x 10 s 



(Model HiCR), close to the value inferred by lWebberl (jl998[ ) for the local interstellar medium 



In addition, we consider the effect of varying the core rotation rate Qq- Our reference 
value Qq = 10~^^ corresponds to a ratio of rotation al to gravitational bin ding energy 



of /3 = 0.025, which is typical of the cores discussed in iGoodman et al.l ( 119931 ). A spread 
exists for the inferred /3 values, which motivates us to consider two additional rotation rates: 
fio = 5x10-1^ (Model LoROT) and 2x10-^^ g-i (Model HiROT), corresponding to /3 = 0.006 
and 0.1, respectively. 

The different variants of the reference model are compared in Fig. [TJ which plots the 
infall and rotation speeds on the equator at a common time t = 5.5 x 10^^ s for all models 
except Model HiROT (which has yet to form a central object at this time, because the 
collapse is significantly delayed by the combination of a moderately strong magnetic field 
and fast rotation). The most striking feature of Fig. [7] is that the rotation speed is essentially 
zero at small radii (within ~ 10^^ cm) for all of the models shown. The same is true for Model 
HiROT (not shown in the figure) at later times (greater than ~ 6.2 x 10^^ s"^); when the fast 
rotating core has collapsed. Apparently the magnetic braking is strong enough to remove 
essentially all of the angular momentum of the equatorial material inside the AD shock for 
the realistic ranges of grain size distribution, cosmic ray ionization rate and core rotation 
rate explored here. 

Suppression of RSD in some cases is easy to understand. The elimination of small 
grains in Model LG and the higher cosmic ray ionization rate in Model HiCR increase the 
densities of electrons and ions, which strengthens the coupling between the magnetic field 
and the bulk neutral matter. Since a better magnetic coupling is expected to make the 
braking more efficient, it is not surprising to find that the RSD remains suppressed as in 
the reference model. In these two cases, the AD shock is located at a slightly smaller radius 
compared with the reference case (see Fig. [7]). We interpret this effect as a result of the 
better coupling, which forces the bulk neutral material to collapse more slowly, which in 
turn leads to a somewhat lower central mass (8.62 x 10^^ g for Model LG and 8.69 x 10^^ g 
for Model HiCR) than in the reference model at the same time (9.29 x 10^^ g). Associated 
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Model REF 




^Q14 ^Q15 ^Q16 ^Q17 

Radius along equator (cm) 

Fig. 7. — Comparison of the infall and rotation speeds along the equator for variants of 
the reference model with different grain size distribution (Model LG), cosmic ray ionization 
rate (Models LoCR and HiCR), and rotation rate (Model LoROT), at a representative time 
t = 5.5 X 10^^ s. The nearly zero rotating speed at small radii indicates the absence of a 
rotationally supported disk due to efficient magnetic braking for all cases shown. 
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with the lower mass is a lower magnetic flux accumulating outside the central object, which 
drives a smaller AD shock. 

The situation with Model LoCR is the opposite. The lower cosmic ray ionization rate 
decreases the degree of ionization in the core, which weakens the coupling between the 
magnetic field and the bulk neutral matter. The weakening of magnetic coupling is expected 
to reduce the braking efficiency in principle. The braking is not weakened enough, however, 
to enable an RSD to form, as evidenced by the vanishing rotation speed and fast infall inside 
the AD shock. In this case, the shock radius is somewhat larger than that of the reference 
model (see Fig. [7]), because of a faster post-shock infall, which leads to a larger central mass 
(9.76 X 10^^ g for Model LoCR) and thus a larger left-behind magnetic flux for shock driving. 

Neither Model LoROT nor HiROT produced a rotationally supported disk. The former 
is to be expected, since the angular momentum is lower for a more slowly rotating core, 
and should be easier to remove by magnetic braking. The slower rotation also presents a 
lesser obstacle to the collapse, allowing the formation of a larger central mass (1.24 x 10^^ g) 
and a larger AD shock (see Fig. [7]) than in the reference model at the same time. In the 
opposite case of HiROT, the stronger rotational support delayed the formation of a central 
object until a much later time. Once the central object has formed, magnetic braking is 
again strong enough to remove essentially all of the angular momentum from the material 
near the central object, preventing an RSD from forming during the protostellar accretion 
phase. 

We conclude that disk formation is suppressed by a moderately strong magnetic field in 
the presence of ambipolar diffusion for the range of parameters that we consider realistic. 

3.3. Ohmic Dissipation 

How does the Ohmic term in the induction equation (1) affect the core collapse and 
disk formation? To address this question, we repeat the reference model, but with the 
Ohmic dissipation included (Model REF^o in Table 1) in addition to the ambipolar diffusion 
(although not yet the Hall term, see below). It turns out that the Ohmic dissipation modifies 
the dynamics of the core collapse and protostellar accretion relatively little, as can be seen 
from Fig. [HI which compares the equatorial infall and rotation speeds of Model REF^o and 
the reference model at four representative times. The speeds for the two cases are barely 
distinguishable over most of the space. The conclusion is that in the presence of ambipolar 
diffusion, Ohmic dissipation does not enable the formation of RSD. 

The reason that Ohmic dissipation is ineffective in modifying the flow dynamics in 
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^Q14 ^Q15 ^Q16 ^Q17 

Radius along equator (cm) 

Fig. 8. — Comparison of the infall and rotation speeds along the equator for Model REFyio 
(with both ambipolar diffusion and Ohmic dissipation; dotted lines) to the reference model 
(with only ambipolar diffusion; solid), at representative times t = 5 x 10^^, 5.5 x 10^^, 6 x 10^^, 
and 6.5 x 10^^ s. The similarity between the two sets of curves shows that Ohmic dissipation 
does not affect the flow dynamics much. 
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the presence of ambipolar diffusion is simple: the Ohmic diffusivity rjo is smaller than the 
ambipolar diffusivity 77^ by more than an order of magnitude (see Fig. [9]). The relatively 
small Ohmic diffusivity is a result of the relatively low density in the computational domain, 
with p < 10~^^gcm~^ (or uh ^ 4 x 10^° cm~^) typically. The moderate density is a result of 
strong magnetic braking, which yields a relatively low-density, collapsing pseudodisk (rather 
than a denser, well supported RSD) in the equatorial region. 

Even when ambipolar diffusion is turned off, we do not find any rotationally supported 
disk within our computation domain (outside a radius of 10^^ cm) that is resulted from 
Ohmic dissipation (Model REFo in Table 1). Indeed, the overall flow dynamics in the 
pure Ohmic case is strikingly similar to the reference case that has only ambipolar diffusion 
(compare Fig. dO] to Fig. H]). The reason for the similarity is that Ohmic resistivity does 
not destroy the net magnetic flux that passes from one hemisphere to the other through the 
equatorial plane (because the poloidal field lines are tied to the low density regions well above 
an d below the equ atorial plane where Ohmic dissipation is negligibly small; see discussion 



m 



Shu et al.ll2006l ). Rather, it enables the poloidal field lines to diffuse radially outward, 
similar to ambipolar diffusion. As more and more matter accretes across the field lines into 
the center, the left-behind magnetic flux piles up at small radii, as shown in the left panel of 
Fig. [ini The flux pile-up is qualitatively similar to the pure- AD reference case (see the left 
panel of Fig. H]). It leads to a strongly magnetized equatorial region where the infall speed 
decreases to well below the free-fall value and the rotation is almost completely braked (see 
the right panel of Fig. [TUl) . We have repeated the Ohmic dissipation-only calculation with 
a reduced radius for the inner boundary (from 10^^ cm to 3 x 10^'^ cm or 2AU), and found 
the same result: namely, a realistic level of (classical) Ohmic resistivity is not large enough 
to enable the formation of rotationally suppo rted disks larger than sever al AUs, at least in 



2D (axisymmetry). This is in agreement with iKrasnopolsky et al.l (120101 ). who showed that 



enhanced resistivity is needed for such disks to form. 



3.4. Hall Effect 

Under the conditions encountered in our core collapse calculations, the magnitude of the 
Hall diffusivity in the induction equation (1) is typically larger than the Ohmic diffusivity (see 
Fig. [3]). We therefore expect the Hall term to have a larger effect on the collapse dynamics 
than the Ohmic term discussed in the last subsection. We find that this is indeed the case. 
Fig.[TT]plots the infall and rotation speeds on the equator for Model KEF^^jq, which includes 
all three nonideal MHD terms in equation (1) and an initial magnetic field that points in a 
direction opposite to the initial rotation axis (in the negative "z" direction, and hence the 
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Fig. 9. — Comparison of the Ohmic diffusivity (dashed hne) and the ambipolar diffusivity 
(sohd) along the equator for Model REF^o (with both ambipolar diffusion and Ohmic dis- 
sipation) at a representative time 6 x 10^^ s. The large difference between the two indicates 
that the Ohmic dissipation is relatively unimportant. 
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Fig. 10. — Density map with field lines (left panel) and equatorial infall and rotation speed 
(right) for Model REFo (with Ohmic dissipation only) at t = 6 x 10^^ s. Also plotted for 
comparison in the right panel are the Keplerian speed based on the central mass (upper 
dashed curve), free fall speed (lower dash-dotted curve), sound speed (horizontal dashed 
line), and zero speed line. Note the rapid deceleration of infalling material near the radius 
~ 2 X 10^^ cm. It is caused by magnetic flux pileup at small radii due to Ohmic diffusion, 
similar to the AD shock shown in Figs. H] and [51 
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superscript "-" in the model name). Plotted are the speeds during the transition from the 
pre-stellar core evolution to the protostellar accretion phase, as in Fig. [6] for the reference 
(AD only) model. Comparing Figs. [TT] and M reveals that the Hall term has relatively little 
effect on the equatorial infall speed. In both cases, as more and more mass accumulates 
at the center, an AD shock develops where the collapsing material slows down temporarily, 
before reaccelerating towards the origin. The effect on the rotation speed is much more 
pronounced. The Hall effect enabled the equatorial material in the post-AD shock region to 
rotate faster compared to the reference case. The difference is especially clear at later times, 
when the rotation inside the AD shock is almost completely braked (or even reversed) by 
magnetic braking in the reference case. In the presence of the Hall effect, the post-shock 
material rotates at a speed as high as 5 times the sound speed. The rotation speed remains 
much smaller than the Keplerian speed, however, indicating that a rotationally supported 
disk is not enabled by the Hall effect. The lack of an RSD is also evident from the fact that 
the post-shock region collapses at a high speed at small radii. 

The significant Hall effect on the post-shock rotation speed can be understood as follows. 
The Hall effect depends on the current density j (oc VxB; see eq. [1]), which is predominantly 
in the toroidal direction because of strong magnetic field pinching in the equatorial region. 
It is particularly strong in the equatorial post-shock region, where the field is strong and is 
bent outward significantly (see the left panel of FigJl]). The toroidal current drives a twist 
of the field lines in the azimuthal direction, which in turn generates a torque that acts to 
spin up the post-shock material. 



An interesting feature, pointed out for example in the disk-wind study of lWardle fc Konigl 



(119931 ) ■ is that when the field direction is reversed , the torque induced by the Hall effect 



changes direction as well (IKrasnopolsky et al.ll201ll ). This is illustrated in Fig. [121 Whereas 



the reference model has a post-shock region nearly completely braked in the absence of the 
Hall effect, both Model REFJ^^^q and Model REF^^^q have substantial, supersonic rotation 
in the post-shock region, although in opposite directions. In particular, in Model KEF^^^q 
where the initial magnetic field is aligned (rather than anti-aligned) with the initial rotation 
axis, the magnetic torque induced by the Hall effect has forced the equatorial material at 
small radii to rotate in a direction opposite to the material at larger distances; the resultant 
shear may induce instabilities in 3D that should be investigated in the future. The change 
in the direction of the torque from Model REF^^^q to Model REF^ji^q is due to the change 
in the poloidal field direction which, for the same outward bending of field lines, produces 
a flip in the direction of the toroidal current. Nevertheless, in neither field orientation was 
the Hall spin-up strong enough to produce a rotationally supported disk, as evidenced by 
the rapid infall speed at small radii where the rotation speed is relatively high (but still well 
below the local Keplerian speed). Indeed, there is relatively little change in the infall speed 




Fig. 11. — Infall and rotation speeds along the equator for Model REF^^jq (including all 
three non-ideal MHD terms) during the transition from the pre-stellar core evolution phase 
to the protostellar accretion phase, between t = 4.4 x 10^^ and 4.6 x 10^^ s, as in Fig. |6l The 
Hall effect enables the material at small radii inside the AD shock to rotate at supersonic 
values. 
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with or without the Hall effect, because the radial current density is relatively small, making 
the Hall effect less important in that direction. 

We will return to the Hall effect in section El where we consider the simpler case of the 
collapse and spin- up of an initially non-rotating core due to Hall effect. 

4. Weak Magnetic Fields and Potential Disk Formation 

We now consider cases with a weak initial magnetic field of Bq = 10.6 /iG (at an initial 
molecular hydrogen density of 71^2 = lO^cm"^), which is 30% of the reference value. It 
corresponds to a dimensionless mass-to-fiux ratio of A = 9.73 for the core as a whole, and 
Ac = 14.6 on the central magnetic flux tube. If the dense core were to condense more or less 
isotropically out of a more diffuse material of lO'^'^ cm~^ in density under the flux freezing 
condition, it would require only an unrealistically weak field of 1.06 fiG in the diffuse gas for 
the core to have Bq = 10.6 /^G. For this reason, we believe that this value of Bq is probably 
as low as, if not lower than, the minimum field strength that can be reasonably expected in 
the type of dense cores under consideration. 

Weaker fields are expected to be less efficient in magnetic braking. This is because 
the braking rate involves the product of the toroidal and poloidal field strengths, and thus 
generally scales with the field strength as B^. Indeed, for Model WREF which has the same 
parameters as the reference model (Model REF) except for a weaker field of Bo = 10.6 /xG, a 
small (~ 20 AU) rotationally supported disk is formed temporarily early in the protostellar 
accretion phase, around the time t ~ 3.68 x 10^^ s, when the central mass is only ~ 0.07 Mq. 
This disk disappears at later times, however, because of strong magnetic braking, which 
drives a powerful outfiow before the disk disappears (see Shang et al., in preparation). 

The suppression of the rotationally supported disk at later times is illustrated in Fig. 
[T51 It includes a snapshot of the collapsing core (left panel) and a plot of the equatorial 
infall and rotation speeds (right panel), at t = 5 x 10^^ s when the central mass is 0.29 Mq. 
At this time, the prominent polar outfiow at the earlier times has disappeared. It is replaced 
by a polar region of strong infall, mostly along the magnetic field lines. The lack of outfiow 
indicates the rotating disk that drives the outfiow no longer exists. The disk suppression is 
shown clearly in the equatorial rotation speed, which is close to zero inside the AD shock 
(around 10^^ cm). The rapid, supersonic infall at small radii further supports the lack of an 
RSD, as in the more strongly magnetized standard model (see Fig. |l]and|5]). 

Disk suppression is not unique to Model WREF. In the left panel of Fig. [H] we show 
three additional examples (Models WLG, WHiCR and WLoROT in Table 1), where the disk 
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Fig. 12. — Infall and rotation speeds along the equator for two models of opposite ini- 
tial magnetic orientation, Model REF^^jq (dashed lines) and Model REF^^^q (dotted), at 
t = 4.55 X 10^^ s. The reference model (solid) without the Hall effect is also plotted for 
comparison. 
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Fig. 13. — Density color map, velocity field (white arrows) and magnetic field lines (yellow) of 
the collapsing core (left panel) and equatorial infall and rotation speeds (right panel) for the 
standard weak magnetic field model (Model WREF) at a relatively late time t = 5 x 10^^ s, 
showing that the RSD formed earlier has by now completely disappeared. Also plotted for 
comparison in the right panel are the Keplerian speed (upper dashed) and free fall speed 
(lower dash-dotted) based on the central mass, the sound speed (horizontal dashed hne), 
and zero speed line. The "effective ion speeds" defined in eq. (4) are shown as dotted lines. 
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formation is suppressed at a relatively late time t = 5 x 10^^ s, when the central mass is 0.40, 
0.34 and 0.65 Mq, respectively. These models have the same weak magnetic field as Model 
WREF, but have either a larger grain size (Model WLG), a higher cosmic ray ionization rate 
(Model WHiCR), or a lower initial rotation rate (Model WLoROT). The lack of disk is not 
too surprising for Models WLG and WHiCR, because the magnetic fields are better coupled 
to the neutral matter in these models than in Model WREF. The absence of a disk in Model 
WLoROT is also expected, because its more slowly rotating core is more easily braked. 

Some rotationally supported disks do form with other choices of parameters, however, 
at least at early times. The disk formation is illustrated in the right panel of Fig. [TU where 
the equatorial infall and rotation speeds for Models WLoCR (same as the standard weak 
field model but with a lower cosmic ray ionization rate), WHiROT (with a higher initial core 
rotation rate) and VWREF (with an unrealistically low Bq = 3.54 fiG) are plotted at a early 
time, when the central mass is only 1.79 x 10^^, 1.27 x 10^^ and 1.05 x 10^^ M©, respectively. 
There is more mass in the rotationally supported disk than at the center, which is why the 
rotation curve is non-Keplerian except close to the origin. Each of the disks drives a strong, 
sometimes chaotic, outfiow, which makes it hard to continue the non-ideal MHD simulation 
reliably to much later times. The disk in the weakest field case (Model VWREF) evolves 
into a ring, which may fragment in 3D. It is unclear whether the early disk in the other two 
cases can survive subsequent magnetic braking or not (it did not in the standard weak field 
model WREF). Once a disk becomes self-gravitating, gravitational torque will likely become 
important in the disk dynamics. This important effect is not captured in our axisymmetric 
simulations. In any case, it is clear that disk formation in the moderately weak magnetic field 
case is more complicated than the moderately strong field case, with the outcome depending 
on the core rotation rate and the degree of field-matter coupling, at least at early times. 
Paradoxically, the weaker field cases are more difficult to simulate because of strong, chaotic 
outfiows. More work is needed before firmer conclusions on disk formation can be drawn. 



5. Hall Spin- up of Non- Rotating Envelope 



Of the three non-ideal MHD effects, the Hall effect is the least explore d in the context 
of co re collapse and disk formation (see, however, iKrasnopolsky et al.ll201ll and lC. Braiding 



201ll ). We have seen in §3.41 that it can spin up the nearly completely braked post-AD shock 



material in an initially rotating core to a significant speed, and the sense of the Hall-induced 
rotation depends on the orientation of the magnetic field. This Hall spin-up can be illustrated 
even more clearly in the collapse of an initially non-rotating core, where any rotation that 
develops subsequently must come solely from the Hall effect. 
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Fig. 14. — Left panel: Infall and rotation speeds on the equator for those models with a 
weak magnetic field that do not have a rotationally supported disk at a relative late time 
t = 5 X 10^^ s. Right panel: Same as the left panel but for those models that do have what 
appear to be the rotationally supported disks at early times, when the central mass is only 
~ 1% of the core mass. Whether these (small) disks can survive to much later times or not 
is unclear. 
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We will concentrate on Model NoROT^^^ where all three non-ideal MHD effects are 
included and the initial magnetic field is anti-parallel to the rotation axis. We have confirmed 
that the case with opposite field orientation (Model NoROT^^^) produces identical results, 
except that the sign of the Hall-induced rotation is flipped, as expected. The left panel of 
Fig. [15] gives an overall impression of the Hall-induced rotation on an intermediate scale 
of ~ 10^^ cm, at a representative time t = 4.4 x 10^^ s, when the central mass is 0.26 Mq. 
Note the alternating pattern of negative and positive rotation speeds, with maximum values 
reaching ~ lO^cms"^, which is much higher than the sound speed (2 x lO^cms"^). On 
smaller scales, the structure is dominated by a dense flattened pseudodisk that is visible 
from the iso-density contours in Fig. [151 

The dense equatorial pseudodisk is collapsing as well as spinning in the positive az- 
imuthal direction, as shown pictorially in the right panel of Fig. [TJ] and more quantitatively 
in the left panel of Fig. [161 The collapse shows rapid deceleration near ~ 2 x 10^^ cm and 
reacceleration interior to it, characteristic of an AD shock. The infall speed is similar to 
that in the pure- AD case, indicating that the Hall effect modifies relatively little the overall 
collapse dynamics. The difference in rotation speed is more pronounced, especially in the 
post-shock region, where the maximum rotation speed exceeds 10^ cms~^, comparable to the 
peak value on the large scale shown in Fig. [151 The rotational component of the "effective 
ion speed" (defined in equation [4]) is larger than the neutral rotation speed, indicating that 
a magnetic force is exerted in the positive azimuthal direction to drive the ion- neutral drift. 
It is the same force that torques up the pseudodisk. The force is particularly large in the 
postshock region, because the magnetic field is strong there. Nevertheless, the spin up fails 
far short of reaching a Keplerian rate, which is the reason why the overall collapse dynamics 
is little affected. 

In an initially non-rotating core, any spin up in one direction must be offset by a spin 
up in the opposite direction, so that the total angular momentum is conserved. In the region 
that extends from the origin up to the AD shock, the net angular momentum is positive (see 
the right panel of Fig. [T6|) . because it is dominated by the positively spinning pseudodisk. 
Outside the AD shock, the net angular momentum is negative, dominated by the hour-glass 
shaped counter-rotating region shown in left panel of Fig. [TJl One may expect the positive 
and negative angular momenta to sum up to zero over the entire computational volume. 
However, this is not the case, because the total angular momentum is dominated by the 
material at large distances (near the core edge) whose rotation speed is small but non-zero 
(see the left panel of Fig. [TB]) . Because the mass at larger distances is larger and has a longer 
lever arm, it dominates the total angular momentum. The total positive angular momentum 
inside the computation domain is 6.52 x lO^^gcm^s"^ whereas the total negative angular 
momentum is —9.38 x lO^^gcm^s"^. They do not cancel out exactly. Some of the angular 
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Fig. 15. — Left panel: map of the Hall-induced rotation speed during the collapse of an 
initially non-rotating core (Model KEF^^jq) at a representative time t = 4.4 x 10^^ s, when 
the central mass is 0.26 M0. Also plotted are magnetic field lines (yellow), isodensity contours 
(red, spaced by a factor of 10), and velocity vectors (white). Right panel: map of log{p) 
(color) and velocity field (white vectors) in the equatorial plane, showing the spin up of the 
dense, post-AD shock region due to Hall effect. 
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Fig. 16. — Left Panel: equatorial infall and rotation speeds of the neutral matter (solid lines) 
for the collapse of an initially non-rotating core. Also plotted for comparison are "effective 
ion speeds" (defined in equation [4], dotted lines), and the infall and rotation speeds in an 
AD-only case that does not have any Hall effect (dashed lines). Right panel: Positive (solid 
line) and negative (dashed) angular momentum per unit radius dL/dr as a function of radius. 
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momentum must have left the computation box, through torsional Alfven waves. Dividing 
the total net angular momentum by the total mass left in the simulation box (1.47 x 10'^'^ g) 
yields 3.80 x 10^^ cm^ s~^, which is the average specific angular momentum. For a core size of 
10^^ cm, the corresponding characteristic rotation speed is 3.8 x 10^ cms~^, about 2% of the 
sound speed, much smaller than the rotation speed achieved on the 10^^ cm scale or smaller. 
We conclude that, despite the localized, supersonic rotation induced on small scales, the 
influence of the Hall effect on the global dynamics is limited. 



6. Discussion 

6.1. Why is Protostellar Disk Formation Difficult in Magnetized Cores? 

6.1.1. Ideal MHD Limit 

The fundamental reason for the difficulty in forming protostellar disks in magnetized 
cores is that the protostellar collapse concentrates magnetic flux at small radii, precisely 
where the rotationally supported disk (RSD) tends to form in the absence of magnetic 
braking. The basic difficulty can be seen most clearly in the ideal MHD limit, where magnetic 
flux is dragged into the central star (because of flux freezing) to form a split magnetic 
monopole. The rapid increase in fleld strength towards the central object enables the (split) 
monopolar field t o brake the circumstellar disk catastrophically, as shown analytically by 
Galli et al.l ( 120061 : see their Fig. 1 for a sketch of the expected field geometry). 



A case can also b e made for catastrophic disk braking in the ideal MHD limit from 
numerical simulations (lAUen et al.l l2003t iMellon fc Lill2008l : iHennebelle fc Fromang||2008l ). 
although it is less clean cut. The reason is that, as the mass of the protostar grows, more 
and more magnetic fiux is dragged to the origin, which creates a stronger and stronger 
(split) magnetic monopole that squeezes more and more strongly on the material on the 
equator from above and below. When the oppositely directed magnetic fields above and 
below the equator are squeezed within a single cell of each other, numerical reconnection 
becomes unavoi dable. The expec ted numerical reconnection is present in the ideal MHD 



simulations of iMellon fc Lil ( 120081 ). especially for relatively strongly magnetized cores (with 
a dimensionless mass-to-fiux A of several or less; see their Figs. 14-16), where no RSDs form 
but the numerical results are complicated by reconnections. The collapse of more weakly 
magnetized cores of A ~ 10 is not significantly affected by reconnections, and yet no RSDs 
form either. The conclusion from the ideal MHD simulations is that magnetic braking is 
efficient enough to suppress disk formation for A < 10. 
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The above conclusion is strengthened by the resistive MHD disk-formation calculations 



of iKrasnopolsky et al.l (120 lOf ). which included a wide range of (prescribed) resistivity. They 
found that, as the resistivity is decreased below a certain value, (numerical) reconnection 
starts to become important (as expected), which complicates the interpretation of the nu- 
merical results. However, before the reconnection sets in, the RSD is already completely 
braked by a moderately strong magnetic field. Extrapolating the results for those clean, low 
resistivity runs without numerical reconnection to the zero resistivity limit (where numerical 
reconnection is unavoidable) indicates that complete disk suppression also holds true for the 
ideal MHD case. 



Machida et al.l (120 lOf ) appears to have come to a different conclusion. They found a 
10^-AU scale disk in their 3D nested-grid simulations in the ideal MHD limit (their Model 
4), even though their core is strongly magnetized (with a global dimensionless mass-to-fiux 
ratio of A = 1, see their Table 1). It is unclear whether the disk is rotationally supported or 
not. If yes, the result would be hard to understand. For such a strongly magnetized core, 
the mass accumulation at the protostar should produce a strong (split) magnetic monopole 
in the ideal MHD limit, whic h is expected to tr igger powerful numerical reconnection, as 



discussed above and shown in lMellon fc Lil (120081 ): we have re-run their Model 4 and found 



the expected episodic reconnections and no rotationally supported disk. There were no 
reconnections mentioned in their paper, and the apparent lack of reconnection may be an 
indication that considerable numerical diffusivity acts to reconnect magnetic field lines of 
opposite polarities efficiently and prevent magnetic flux from accumulating near the protostar 
to form the expected split monopole in the first place. Since the trapping of magnetic flux 
at small radii lies at the heart of the efficient braking that renders disk formation difficult 
in the ideal MHD limit, inability to do so numerically may weaken the braking efficiency 
artificially and lead to disk formation. 



6.1.2. Non-Ideal MHD Effects 

In lightly ionized dense cores of molecular clouds, non-ideal MHD effects are to be ex- 
pected. Non-ideal effects, particularly ambipolar diffusion and Ohmic dissipation, enable the 
bulk neutral matter to move across magnetic field lines, breaking the flux freezing condition 
that is responsible for the formation of the central split magnetic monopole which, in turn, 
is responsible for the catastrophic disk braking in the ideal MHD limit. The elimination of 
the central split monopole does not necessarily mean, however, that the magnetic braking 
would automatically be weakened enough for an RSD to appear. The reason is that the 
magnetic flux that would be trapped in the central split monopole in the ideal MHD limit is 
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now concentrated in a small, but finite, circumstellar region ii istead, as first demonstrated by 
Li fc McKed fll996[ ) in the case of ambipolar diffusion (AD). iKrasnopolsky fc Konigll fj2002l ) 
showed semi-analytically that the AD-induced flux concentration at sm all radii can in prin- 
ciple suppress disk formation completely, just as in the ideal MHD limit. iMellon fc Lil (j2009[ ) 
showed numerically that RSDs are indeed suppressed by a moderately strong magnetic fleld 
(with A ~ several) in the presence of ambipolar diffusion for a reasonable range of cosmic 
ray ionization rate. 



The result of iMellon fc Lil (120091 ) is strengthened by the simulations presented in §3.11 
We improved over their calculations by self-consistently computing charge densities including 
dust grains and by extending the computation to the prestellar phase of core evolution leading 
up to the central mass formation. The extension allows us to explore the angular momentum 
evolution and disk formation during the transition between the prestellar and protostellar 
phase of star formation. We flnd that the tendency to form a disk is stronger around the 
time of initial protostar formation than at later times (see Fig. [6]). This is because there is 
as yet little magnetic flux accumulated near the center and it takes time for the magnetic 
braking to remove angular momentum. Once enough magnetic flux has accumulated near the 
protostar to drive a well-developed accretion (C-)shock, the fleld strength in the post-shock 
region is typically strong enough to remove most of the angular momentum of the material 
falling into the central object, as long as the dense core is moderately strongly magnetized 
to begin with (with A ~ a few to several). 

The conclusion that the RSD is suppressed by a moderately strong magnetic fleld in 
the presence of AD is robust, because the size of the AD shock and the post-shock fleld 
strength are rather insensitive to cosmic ray ionization rate and the grain size distribution 
(see Fig. [7]). They are determined mostly by the global requirements that (1) most of the 
magnetic flux associated with the central stellar mass be redistributed in the post-shock 
region, and (2) the strong p ostshock magnetic field be confined by the infall ram pressure 
(see equations [8] and [10] of lLi fc McKedll996l ). as long as the width of the C-shock is small 
compared to the radius of the shock. In principle, if the ionization level is decreased by a 
arbitrarily large factor, the magnetic field would eventually decouple completely from the 
bulk neutral material and an RSD would form. In practice, however, the RSD is suppressed 
even for the highly conservative case of both an unrealistically low cosmic ray ionization rate 
of ( = 10^^^ s^^ and an MRN grain size distribution that contains a large amount of small 
grains (Model LoCR in Table 1). Both the grain growth expected in dense cores and a more 
realistic (higher) cosmic ray ionization rate tend to make the magnetic field better coupled 
to the bulk neutral material and the magnetic braking more efficient. 



The Ohmic dissipation does not change the above picture much, because the Ohmic 
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diffusivity rjo is well below the ambipolar diffusivity tjad for the density range uh ^ 10^^ cm ^ 
that is crucial for disk formation. We can estimate the ratio of the two diffusivities through 

Vo 1 



Vad rie /3c/3e 

where the subscripts "c" and "e" denote, respectively, the charged species whose contribution 
dominates the AD term and the electrons that are mainly responsible for the Ohmic term. 
The (dimensionless) Hall parameter P = ru (where r is the coUisional damping time of the 
motion of a charged species relative to the neutral and u the cyclotron frequency) provides 
a measure of how well a charge is tied to the magnetic field. The relative unimportance of 
Ohmic dissipation comes mainly from the fact that electrons are extremely well tied to the 
magnetic field in the density regime of interest, with /3e ^ 1. For example, at a representative 
density of uh = 3 x lO^cm"^, we find /3e ~ 2 x 10^ for the field strength-density relation 
given by equation (5). In the large grain (LG) case that we have considered, the AD term is 
dominated by metal ions, M"*". Their number density is close to the electron number density 
{nM+ ~ ne) and they are well tied to the magnetic field at the representative density, with 
/3m+ ~ 40. From equation ([6]), we expect tio/tiad ~ 10^^, which is in agreement with 
the computed values shown in the right panel of Fig. [31 In the case of MRN grain size 
distribution, the AD term is dominated by the small negatively charged grains, g~ . At the 
representative density, there are about 10^ small negatively charged grains for each electron 
{rig- lO^ne) and the small charged grains are marginally tied to the magnetic field fig- ~ 1. 
Both the smaller electron abundance and the weaker coupling of the grains to the field tend 
to make Ohmic dissipation more important relative to ambipolar diffusion. Nevertheless, the 
electrons are so well tied to the magnetic field that, even in this case, the Ohmic diffusivity is 
still much smaller than the ambipolar diffusivity, by a factor of ~ 2000 at the representative 
density. This estimate is again in agreement with the computed values shown in the left 
panel of Fig. [3l 

The Hall effect is expected to be more important than Ohmic dissipation in diffusing 
the magnetic field in the density regime under consideration. This can be seen from the ratio 
of Hall and ambipolar diffusivities 

^ ~ - (7) 

which is applicable under the conditions that the same charged species "c" dominates both 
the Hall and AD terms and /3c > 1. These conditions are satisfied for the large grain (LG) 
case at the representative density nu = 3 x 10* cm""^, where metal ions dominate both terms 
and ~ 40. In this case, the Hall diffusivity is only about 2% of the ambipolar diffusivity 
(but still much larger than the Ohmic diffusivity). It becomes comparable to the ambipolar 
diffusivity in the case of MRN grain size distribution, where both terms are dominated by 
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small dust grains that are marginally coupled to the magnetic field (with a Hall parameter 
/3 of order unity) at the representative density. At higher densities, the Hall diffusivity 
is expected to exceed the ambipolar diffusivity as small grains become even less well tied 
to the magnetic field. However, the contribution from positively charged grains start to 
cancel out that from negatively charged grains, leaving the Hall diffusivity comparable to 
the ambipolar diffusivity over a wide range of density (see the left panel of Fig. [3]). The Hall 
effect therefore does not increase the magnetic diffusivity by more than a factor of a few. As 
such, it is not expected to greatly change the global flow dynamics, especially the structure 
of the ambipolar diffusion-induced accretion shock, which lies at the heart of the magnetic 
braking catastrophe. It does, however, introduce a new ingredient into the problem: it can 
actively torque up a magnetized col lapsing envelope, eve n if the envelope is non-rotating to 
begin with, as flrst pointe d out b y IWardle fc Na (19991 ) and demonstrated numerically in 
Krasnopolsky et al. ( 2011 : see also C. Braiding 20 111 ). On the scale of several AUs or larger 
that we can resolve in our non-ideal MHD simulations ( §3.41 and |5|). the Hall spin- up does 



not reach the Keplerian speed. The angular momentum gained through Hall spin-up may, 
however, be conducive to the formation of RSDs on smaller scales, particularly at high enough 
densities where electrons begin to decouple from the magnetic field a r id Ohmic dissipatio n 
becomes the dominant process for field diffusion ( iMachida et al.ll2010l : iDapp fc Basijl2010l ). 
Nevertheless, the problem of catastrophic magnetic braking that prevents the formation of 
a sizable RSD of tens of AUs or larger is not resolved through the three non-ideal MHD 
effects. 



6.2. Limitations and Future Directions: How to Form RSDs? 

We are unable to produce robust, large-scale, rotationally supported disks in our non- 
ideal MHD simulations for dense cores magnetized to a realistic level. And yet, rotationally 
supported disks are observed around young stars, at least at relatively late times, after the 
massive envelope has been removed (whi ch reveals the embedded disk for direct rotation 
measurement; see IWilliams fc Ciezall201ll for a review). Clearly, one or more assumptions 
made in our calculations must break down in order for the observed (late-time) RSDs to 
form. These include (1) assumptions made in the setup of the numerical problem because 
of computational constraints and (2) additional physical effects that we have not taken into 
account. We comment on these limitations and their relevance to RSD formation in turn. 
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6.2.1. Numerical Limitations and Possible Ways to Form RSDs 



We have restricted our problem setup to 2D (assuming axisymmetry) , which greatly 
reduced the computational demand. A potential drawback is that the imposed symmetry 
may have enhanced the ability of the AD shock in trapping magnetic flux near the central 
object, which lies at the heart of the catastrophic magnetic braking that prevents disk 
formation. The reason is that the high magnetic pressure may force the trapped fleld lines to 
escape along the path of least resistance when the axisymmetry is broken. It is plausible that 
some magnetic flux loaded with relatively little matter (recall that the post-shock region is 
magnetically subcritical and supported against free-fall collapse by magnetic forces) would 
act as a "light" fluid and escape in some azimuthal directions, allowing less magnetized 
flngers of "heavy" material to sink clos er to the center in other directions. This type of 
interchange instability was considered in iLi fc McKed ( 119961 ). It may weaken the magnetic 
fleld strength at small radii enough to enable disk formation. Investigation of this possibility 
is now underway. 

Another limitation of the current problem setup is that we adopted a relatively large 
central hole around the protostar (typically 10^^ cm or 6.7 AU in radius). We have experi- 
mented with smaller holes by a factor of 2-3 in a few cases and found quantitatively similar 
results. However, it is difficult to reduce the hole size by a much larger factor, because of the 
constraints on the time step, which decreases as the square of the hole size (assuming the 
same resolution in the 6'-direction) in our explicit treatment of the non-ideal MHD terms. A 
drawback is that we are unable to determine whether small, (sub-)AU-scale RSDs can form 
during the main protostellar accretion phase or not. The existence of such small RSDs is 
suggested by the powerful m olecular outflows ubiq uitously observed around deeply embed- 
ded. Class protostars (e.g., iBontemps et al.lll996[ ): they are generally thought to be driven 



by a fast primary wind launche d magnetocentrifugally from the inner part of a Keplerian 
disk close to the central object (IShu et al.ll2000l : iKonigl fc Pudritzl 120001 ). 



Small RSDs can in principle form at high enough densities where electrons start to 
decouple from the magnetic fleld. Before thermal ionization of alkali metals becomes im- 
portant, the Ohmic dissipation can reduce the local current density, making it hard for the 
magnetic fleld to bend in the poloidal plane (which limits the flux accumulatio n at small 



radii) and to twist in the azimuthal direction (which weakens magnetic braking; IShu et al. 



20061 : iMachida et al.l 120071 : iKrasnopolskv et al.ll2010l : iDapp fc Basull2010f ). A worry is that, 
in the presence of only ambipolar diffusion and classical Ohmic dissipation, our (2D ax- 
isymmetric) calculations showed that little angular momentum is left at the inner edge of 
the computation domain (6.7 AU), making the formation of RSDs interior to it difficult. 
However, the Hall effect can spin up the flow that collapses through the inner boundary 
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to a supersonic (although still locally sub-Keplerian) rotation speed, perhaps making the 
formation of small RSDs possible. Alternatively, the 3D magnetic interchange instability 
discussed earlier may weaken the magnetic braking interior to the AD shock enough to facil- 
itate the formation of RSDs in general, and the small RSDs required for fast wind launching 
in particular. When and how the small RSDs grow to large RSDs observed around relatively 
evolved YSOs remains uncertain, and may require additional physical effects that have not 
been investigated in detail in this context to date. 



6.2.2. Additional Physical Effects for RSD Formation 



An important ingredient for low-mass star formation is protostellar o utflow . It may 
play a crucial role in disk formation. As first proposed by iMellon fc Lil ( l2008l : see their 
§6.2.2), the outflow can strip away the slowly rotating protostellar envelope, which brakes the 
equatorial infall material that is magnetically linked to it and that tries to spin up and form a 
rotationally sup ported disk. Part of th e envelope may be removed by the core collapse process 
itself 



isee, e.g.. 



Machida et al.l |2010[) but, if the e fficiency of star fo rmation in individual 



low- mass cores is of order 1/3 (lAlves et al.l 120071 : lAndre et al.ll2010[ ). the majority of the 
envelope mass must be re moved by some other process, most likely a (fast) protostellar wind 
(IMatzner fc McKeell2000l ). Indeed, the bipol ar molecular out f low, thought to be primarily the 
envelope material swept up by a fast wind fIShu et al.lll99ll : IShang et al.ll2006l ). is observed 
to have a narrow jet-like appearance along the axis du ring the early Class phase and the 
opening angle at the base increases as the YSO ages (lArce fc SargentI l2006l ). As the fast 
wind sweeps out an increasingly wider polar region in the envelope, the braking efficiency of 
the remaining equatorial infall region should decrease, perhaps to a low enough value that a 
large-scale rotationally supported disk can form. 

A specific outflow-enabled large-scale disk formation scenario is as follow. We envision 
the early formation of a small (perhaps AU-scale) rotationally supported disk (unresolved 
by the current generation of instruments) during the Class phase, through the processes 
discussed in §6.2.11 Although the small disk can grow gradually through internal angular 
momentum redistribution (perhaps gravitational torques rather than magnetic stresses since 
magnetic decoupliri g is required for the disk to form in the ffist place, see discussion in 
Dapp fc Basull2010l ). we envision rapid growth in disk size (to, say, 100 AU or more) only 
during the late phase of envelope removal, when the braking of the equatorial infall material 
by the envelope is rendered inefficient by outflow stripping. This envelope-depletion induced 
rapid disk growth may occur towards the end of the main protostellar mass accretion phase, 
perhaps during the transition from the Class to Class I phases of (low-mass) star formation. 
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Detailed calculations and high resolution observations, perhaps using ALMA, are needed to 
test this scenario of late-time formation of a large-scale disk. 

Another possibility for large-scale RSD formation is through enhanced magnetic diffu- 
sivity. If the diffusivity is greatly enhanced over the classical microscopic values considered 
in this paper by so me processes, th e magnetic braking may be weakened enough to allow 



i ty ma y enable RSD formation, and this was demonstrated explicitly in 



(120101 ) . The enhancement in resistivity may come from turbulence (IKowal et al.l l2009l ) 



KrasnoDolskv et al. 


Kowal et al. 


2009 





which is observed in dense cores from nonthermal line width, or current- driven instabilities 



flNorman fc Hevvaertslll985[). although these effects are hard to quantify at the present time 
(see, however, ISantos-Lima et al.l l2010l w 'ho have started to quanti fy the so-called "turbu- 
lent reconnection diffusivity"). Similarly, iKrasnopolsky et al.l ( 120111 ) showed that if the Hall 
diffusivity is large enough, it can enable RSDs to form even in initially non-rotating dense 
cores. 

Numerical diffusion may mimic to some ext ent enhanced rnagneti c diffusion of physical 
origins and lead to large-scale RSD formation. iMachida et al.l ( l2010l ) was able to produce 
large-scale disks in strongly magnetized cloud cores in both the ideal MHD limit and with 
a classical (microscopic) value of Ohmic dissipation. However, as we have argued in §6.1.H 
the lack of episodic reconnections expected in the ideal MHD limit indicates a considerable 
numerical diffusion in their calculations. Another indication is that, in the presence of only 
the classical Ohmic dissipation, we do not find any large-scale RSD (Model REFo in Table 
1), for a good reason: the Ohmic diffusivity enables the magnetic flux that would have gone 
into the central object in the ideal MHD limit to accumulate in a small circumstellar region 
where the magnetic braking is particularly efficie nt fsee Fig. [TO|l. as in the AD case. Such 
a magnetic flux accumulation was not obvious in iMachida et al.r s simulations, which may 
again indicate an enhanced magnetic diffusion, either of numerical origin or through 3D 
effects that are not captured by our 2D calculations. 



6.2.3. Weak Core Magnetization and RSD Formation 

Here we comment on the possibility that dense star forming cores may be magnetized 
to different levels, and disks form preferentially in those that are weakly magnetized. Our 
calculations indicate that, in the presence of ambipolar diffusion, the core mass-to- flux ratio 
need to be greater than at least ~ 10 in order for a rotationally supported disk to form 
and survive to late t i mes. Observationally, a mean value of Xios ~ 4.8 ± 0.4 is inferred by 



Troland &: Crutcherl (120081 ) from the line-of-sight field strength for a sample of dark cloud 
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cores. Applying; geometric corrections would reduce the value statistically by a factor of 2-3 
( IShu et al.l Il999[ ). making it unlikely for the majority of dense cores to be magnetized as 
weakly as A > 10. Since the majority of, if not all, young stars (formed out of all dense 
cores) are thought to have an RSD at some point, we consider it unlikely that weak core 
magnetization is the main reason for RSD formation. 



6.3. Observational Implications: Disk vs Pseudodisk 



We should emphasize that, in our simulations, even though large rotationally supported 
disks are difficult to form, highly flattened dense "disk-like" structures are prevalent (see 
Figs, m [TOl [13], and [15]). This is not surprising because, just like rotation, the magnetic field 
can provide anisotropic su pport to the cloud core, allowing matter to settle along field fines 
into fiattened structures (IGalli fc Shul Il993l ) . The fact that there are two types of forces 
in nature that can retard (anisotropically) the gravitational collapse naturally leads to two 
types of fiattened structures: rotationally supported disks and magnetically induced (pseudo- 
) disks. If the dense cores are as strongly magnetized as indicated by the currently available 
observations (with a dimensionless mass-to-fiux ratio of several or smaller, see §1]), then there 
is typically more magnetic energy than rotational energy, and a magnetically induced pseudo- 
disk is just as, if not more, likely to form around an accreting protostar as a rotationally 
supported disk, especially in view of the fact that magnetic braking hinders the formation 
of rotationally supported disks but not pseudodisks. It is therefore premature to conclude 
that dense fiattened structures obs erved around deeply em bedded protostars (such as from 
dust continuum observations, e.g., |j0rgensen et al.l l2009l ) are rotationally supported disks 
rather than magnetically induced pseudodisks; the latter can be just as thin as (perhaps even 
thinner than) rotationally supported disks, because of magnetic compression. To confuse the 
situation further, the pseudodisks can have a substantial rotation as well (just not enough to 
provide the full support against gravity) and may or may not collapse at a high speed (see 
again Figs. [U [TD] [T^ and [T^ . Detailed kinematic information, as well as a knowledge of the 
central mass, are needed to establish whether a dense fiattened circumstellar structure is a 
rotationally supported disk or n ot. High resolution o bservations of the circumstellar magnetic 
field structure, such as those in lGirart et al.l ( 20061). will also go a lo ng way towards testing 
the idea of magnetically induced pseudodisk ( jGoncalves et al.l 120081 ). 
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7. Summary 

We have carried out a set of 2D axisymmetric calculations exploring non-ideal MHD 
effects in magnetic braking and protostellar disk formation in rotating magnetized dense 
cores. Our main conclusions are summarized as follows: 

1. For a realistic magnetic field of moderate strength corresponding to a core mass- 
to-flux ratio A ~ 3-4, the magnetic braking is strong enough to remove essentially all of 
the angular momentum of the material that accretes onto the central object in the presence 
of ambipolar diffusion under a wide range of conditions in 2D. Any large-scale (greater 
than several AUs) rotationally supported disk (RSD) is suppressed by the formation of an 
ambipolar diffusion-induced accretion shock, which traps a strong magnetic field near the 
central object, leading to efficient magnetic braking of the post-shock material. 

2. On scales greater than ~ 10 AU, realistic levels of Ohmic diffusivity do not enable the 
formation of large-scale RSDs, either by itself or in combination with ambipolar diffusion. 
Furthermore, Ohmic dissipation does not necessarily reduce the magnetic braking efficiency. 
It can make the braking more efficient by enabling magnetic flux accumulation at small radii, 
where the field strength is increased, similar to the case of ambipolar diffusion. 

3. The Hall effect can spin up the post- AD shock material to a significant, supersonic 
rotation speed, although the rotation remains too sub-Keplerian to form an RSD for the 
parameter space explored in this work. 

4. For an unusually weak magnetic field corresponding to a core mass-to-fiux ratio 
A > 10, a small RSD often forms early in the protostellar accretion phase, when the central 
mass is still small. In the majority of cases, the RSD disappears at later times, braked 
strongly by the powerful outflow that it drives. In some cases, particularly when the cosmic 
ray ionization rate is unusually low and the core rotation rate is unusually high, the fate of 
the early disk is unknown because the simulation stops early due to numerical difficulty. 

5. We discussed several possible ways to enable the formation of large-scale RSDs: mag- 
netic instabilities in 3D, early formation of small RSDs at high densities, outflow stripping 
of protostellar envelope, enhanced magnetic diffusion and weak core magnetization. The 
more likely of these possibilities are, in our view, the weakening AD shock in 3D through 
interchange instability, which is expected to decrease the field strength (and thus the brak- 
ing efficiency) near the central object, outflow stripping of protostellar envelope, which may 
allow rapid formation of a large-scale RSD during the transition from the deeply embedded 
(Class 0) phase to more revealed (Class I and II) phase of low-mass star formation, and 
enhanced magnetic diffusivity, which may be driven by turbulence-induced reconnections. 
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